Strong gyrotopy in a chiral toroidal medium 
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In this letter, we present the first experimental study of a new chiral metamaterial consisting of 
toroidal wire windings. We show that the metamaterial exhibits three bands of circular dichroism in 
the GHz range. We discuss the response of the structure in terms of multipole moments, including 
the (magnetic) toroidal dipole moment. 



Toroidal, doughnut-shaped structures are ubiquitous 
in nature, appearing on scales which range from the sub- 
atomic [TJ [2] to the astronomical [3 . On the molecular 
level, the torus shape is preferred by numerous biologi- 
cal and chemical macromolecules, such as DNA conden- 
sates [4], proteins [5] and oligosaccharides [6 , to name 
just a few. Toroidal symmetries are also encountered 
frequently in solid-state systems including carbon nan- 
otubes [7] and ferroelectrics [8[ |9] . Although the behav- 
ior of such systems has been studied extensively, their 
interaction with electromagnetic radiation is not well un- 
derstood. Indeed, within the framework of classical elec- 
trodynamics, unusual phenomena associated with viola- 
tion of Lorentz reciprocity [10 and non-radiating config- 
urations [11 have been predicted for toroidal structures 
and their interactions. Nevertheless, such phenomena are 
usually weak p~2j [13] and, therefore, experimental inves- 
tigations are rare. In this letter, we study a new artificial 
chiral medium, originally suggested in [14 , consisting of 
an array of toroidal wire windings in a metamaterial' con- 
figuration. We show that such a metamaterial exhibits 
strong gyrotropic response which is attributed to differ- 
ent terms of its multipole expansion. 

In contrast to artificial gyrotropic media, where hand- 
edness is usually associated with the direction of a "twist 
vector" following a cork-screw law along the helicity axis, 
the situation is more complicated when the structure 
possesses toroidal symmetry. Here, the twist vector ro- 
tates along the torus, and therefore a corresponding di- 
rection can not be defined. However, although no helic- 
ity axis exists, the structure has two well-defined enan- 
tiomeric forms, corresponding to different directions of 
the winding (see Fig. la). Based on this concept, we 
manufactured a chiral toroidal metamaterial with a unit 
cell consisting of four connected square loops that were 
formed by horizontal and vertical copper wire segments 
(see Fig. lb). The resulting windings were embedded in 
dielectric bars with permittivity e = 4.5 — 0.081z. The 
size of the unit cell was 15x15 mm rendering the meta- 
material non-diffracting up to 20 GHz. Transmission 
experiments were performed at normal incidence, from 
2 to 14 GHz, in an anechoic chamber using two broad- 



band horn antennas (Schwarzbeck M. E. model BBHA 
9120D). The transmission spectra were recorded with a 
vector network analyzer. 



(a) 



(b) 




4 ► 

15 mm 



FIG. 1: (a) Left- and right-handed form of a chiral toroidal 
winding. The electromagnetic response of the structure is 
determined by the handedness of the winding and that of the 
incident wave, (b) Top and side view of the fabricated toroidal 
metamaterial unit cell. 

The intensity and phase of the structure's response to 
right and left circularly polarized light are presented in 
Figs. 2a and 2b, respectively. As it can be seen in Fig. 2a, 
experimental results (solid line) are in good agreement 
with finite element numerical simulations (solid circles). 
Two resonant bands of strong circular dichroism can be 
distinguished at around 4.5 and 10 GHz, where the cor- 
responding transmission ratio for orthogonal circularly 
polarized states reaches 5 and 15 dB with bandwidth 
(FWHM) 370 MHz and 1.6 GHz, respectively. More- 
over, while the rejected polarization state experiences 
high losses, the orthogonal polarization state propagates 
through the structure almost unaffected. This is also ev- 
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ident from the behavior of the phase delay presented in 
Fig. 2b. In addition, the rejected polarization state cor- 
responds to a backward wave, while the other to a normal 
propagating wave, since they both have the same sign of 
phase velocity, but different signs of group velocity. This 
can be considered as a plausible indication for the pres- 
ence of negative refraction [15]. Finally, a third band of 
circular dichroism can also be seen at around 7.5 GHz. 
In this case, however, the dichroism is weaker, while both 
polarization states experience significant losses. 
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FIG. 2: (a) Transmitted intensity for right (blue) and left 
(red) circularly polarized light. The circular dichroism is rep- 
resented by the solid black curve (experimental) and the black 
solid circles (numerical). The position of the dichroism res- 
onances is marked by the blue panels, (b) Phase delay for 
right (blue) and left (red) circularly polarized light. Inset: 
Unit cell of the studied metamaterial and direction of wave 
propagation. 

The observed resonances are of geometrical origin and 
their frequency position is controlled by the ratio of the 
total wire length in the metamaterial unit cell over the 
wavelength of the incident wave, while the presence of the 
dielectric has to be taken also into account. In general, 
geometrical resonances will occur whenever, the total 
wire length equals an integer number of half- wavelengths. 
However, due to the symmetry of the structure only an 
even number of nodes is allowed. Therefore, resonances 
will occur whenever the wire length equals an integer 
multiple of the excitation wavelength (even number of 
nodes). At these resonances, the intensity of interaction 
is very different for left and right circular polarization. 
This is illustrated in Figs. 3a & 3b, where energy density 
maps and power flow lines (following the Poynting vec- 
tor) are shown. While one polarization interacts strongly 
with the structure (Fig. 3a) the orthogonal circular po- 
larization propagates almost unaffected. This behavior 
is also demonstrated in Figs. 3b and 3c, where the real 



part of the current intensity along the wire winding is 
plotted for the case of the low frequency (4 GHz) and 
high frequency (10 GHz) resonances. Although the cur- 
rent configuration is very similar for both polarizations, 
resembling standing waves with four (4 GHz) and six 
nodes (10 GHz), its intensity is significantly different 
which leads to the observed dichroism resonances. 
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FIG. 3: Finite element simulations of toroidal metamaterial. 
(a) and (b) show electric density (color maps) and powerflow 
lines at the 4 GHz resonance for right and left circularly polar- 
ized incident wave, respectively, (c) and (d) present induced 
current configuration on the wire winding for right (red) and 
left (blue) circularly polarized light, at the 4 GHz and 10 GHz 
resonances, respectively. 

In the framework of multipole theory, the observed 
gyrotropy arises from the coupling of the electric and 
magnetic dipole moments, but can also include a contri- 
bution from the electric quadrupole moment, since the 
metamaterial is anisotropic [171 HE] • In order to further 
investigate the metamaterial's response, we calculate the 
induced multipole moments [19] under plane wave illumi- 
nation The multipole moments have been calculated by 
the following relations: 



electric dipole : 
magnetic dipole : \ 
electric quadrupole : Qij 



= \j{r x J)idV 
^KriJj + Jir^dV, 



where z, j = x,y,z and compared their strength in terms 
of radiated power: 



el. dipole radiation : 
m. dipole radiation : 
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el. quadrupole radiation : Pq = -£ %°3 \Qij\ 
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as presented in Fig. 4. As expected, the electric dipole 
moment dominates the response of the system at all 
frequencies, while all multipole moments resonate to- 
gether at the 4 GHz and 10 GHz circular dichroism 
bands. More concisely, at the low frequency resonance, 
the magnetic dipole is much stronger than the electric 
quadrupole, indicating that the main contribution to chi- 
rality comes from the electric dipole-magnetic dipole cou- 
pling. This is also supported by the fact that these mo- 
ments are almost collinear as shown in the correspond- 
ing inset of Fig. 4. Similar situation arises also at the 
high frequency dichroism band (10 GHz). However, at 
the weaker 7 GHz resonance, the electric quadrupole 
presents a pronounced resonance, while the magnetic 
dipole moment is weaker. This suggests that the weak 
dichroism at this resonance is a result of electric dipole - 
electric quadrupole coupling. Although ideally the elec- 
tric dipole moment should also be at a minimum, the 
dominant contribution is still coming from the latter, as 
a result of the finite size of the structure. 

In addition, we include in our analysis the (magnetic) 
toroidal dipole moment pQ. Indeed, it has been sug- 
gested that the complete multipole expansion requires 
the inclusion of toroidal moments along with the electric 
and magnetic ones [2T] , While usually these moments 
can be neglected, this is not true for structures compara- 
ble to the wavelength or structures of toroidal symmetry 
[21], as the one considered here. In fact, the toroidal 
moment is known to result in magneto-electric coupling 
[22] , since it enters the multipole expansion in exactly the 
same way as the electric moments [21] . Therefore one can 
expect that coupling between toroidal and magnetic mo- 
ments would lead to optical activity and the question of 
a toroidal contribution to the observed chirality arises. 
In order to quantify the toroidal response of the meta- 
material, we calculate its toroidal moment by the for- 
mula, r = J [(r ■ j)r — r*j]d 3 r and present in Fig. 4 the 
corresponding radiated power P T = " 2 £° 5 \T\ 2 [19 . At 
most frequencies the toroidal moment is comparable with 
the electric quadrupole moment, which suggests that the 
toroidal contribution to the metamaterial response is in 
general no n- negligible. However, as it can be seen from 
the radiated powers, its contribution in the gyrotropy is 
secondary, since the electric and magnetic dipole/electric 
quadrupole moments are much stronger at the dichroism 
resonances. Nevertheless, a much stronger toroidal re- 
sponse is expected at lower frequencies, where electric 
multipole moments vanish and we intend to study this 
regime in the near future. 

In conclusion, we present a first experimental inves- 
tigation of a new type of chiral metamaterials, which 
possesses toroidal symmetry. The studied metamaterial 
shows strong circular dichroism in the GHz frequency 
range as a result of geometrical resonances. The gy- 
rotropic response of the structure can be explained in 
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FIG. 4: Radiated power from the different multipole moments 
of the (numerically simulated) current configuration on the 
wire windings, normalized to the incident power. Insets: Ori- 
entation in the x-y plane of dipole moments at the frequencies 
marked by the corresponding dashed lines. In both cases, the 
z-component of the moments is much smaller than the x and 
y components. 

terms of electric and magnetic dipole as well as electric 
quadrupole moments. Finally, we discussed the toroidal 
response of the system and showed that at specific fre- 
quencies it can be of magnetic dipole/electric quadrupole 
order. 

The authors would like to acknowledge the financial 
support of the Engineering and Physical Sciences Re- 
search Council, UK. 

References 



* Electronic address: N.Papasimakis@soton.ac.uk 
[1] Zel'dovich, Ya. B., Zh. Eksp. Teor. Fiz. 6, 1184 (1958). 
[2] V. M. Dubovik, A. A. Cheshkov, Phys. Part. Nucl. 5, 791 
(1974). 

[3] J. S. Mainstone, Nat. 215, 1048 (1967). 

[4] N. V. Hud and I. D. Vilfan, Annu. Rev. Biophys. Biomol. 
Struct. 34, 295-318 (2005). 

[5] M. M. Hingorani and M. O'Donnell, Nat Rev. Mol. Cell 
Biol. 1, 22-30 (2000). 

[6] W. Saenger, J. Jacob, K. Gessler, T. Steiner, D. Hoff- 
mann, H. Sanbe, K. Koizumi, S. M. Smith, and T. 
Takaha, Chem. Rev. 98, 1787-1802 (1998). 

[7] J. Lium H. Dai, J. H. Hafner, D. T. Colbert, R. E. Smal- 
ley, S. J. Tans, and C. Dekker, Nature 385, 780-781 
(1997). 



4 

[8] A. A. Gorbatsevich and Y. V. Kopaev, Ferroelectrics 

161, 321-324 (1994). 
[9] I. I. Naumov, L. Bellaiche and H. X. Fu, Nature 432, 
737-740 (2004). 
[10] G. N. Afanasiev, J. Phys. D 34, 539-559 (2001). 
[11] G. N. Afanasiev and V. M. Dubovik, Phys. Part. Nucl. 
29, 366 (1998). 

[12] Sawada K and Nagaosa N 2005 Phys. Rev. Lett. 95 
237402 

[13] V. A. Fedotov, K. Marinov, A. D. Boardman and N. I. 

Zheludev, New J. Phys. 9, 95 (2006). 
[14] K. Marinov, A. D. Boardman, V. A. Fedotov, and N. I. 

Zheludev, New J. Phys. 9, 324 (2007). 
[15] J. B. Pendry, Science 306, 1353 (2004). 
[16] R. E. Raab and O. L. de Lange, Multipole Theory in 

Electromagnetism (Oxford, 2005). 
[17] A. D. Buckingham and M. B. Dunn, J. Chemical Soc. A, 

1988 (1971). 

[18] I. P. Theron and J. H. Cloete, IEEE Trans. Antennas 

Propagat. 44, 1451 (1996). 
[19] C. Vrejoiu, J. Phys. A: Math. Gen. 35, 9911 (2002). 
[20] Similar results were obtained under excitation with the 

orthogonal polarization state (y axis of Fig. lb). 
[21] V. M. Dubovik and V. V. Tugushev, Phys. Rep. 187, 145 

(1990). 

[22] M. Fiebig, J. Phys. D 38, R123 (2005). 



